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Direct and large-eddy simulations of the interaction between the wake of a circular 
cylinder and a flat-plate boundary layer are conducted. Two Reynolds numbers are ex- 
amined. The simulations indicate that at the lower Reynolds number the boundary layer 
is buffeted by the unsteady Karman vortex street shed by the cylinder. The fluctuations, 
however, cannot be self-sustained due to the low Reynolds-number, and the flow does 
not reach a turbulent state within the computational domain. In contrast, in the higher 
Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed 
(due, in part, to the higher values of the local Reynolds number Re# achieved in this 
case); some evidence could be observed that a self-sustaining turbulence generation cycle 
was beginning to be established. 


Introduction 

High-lift systems have a significant impact on the 
overall cost and safety of aircraft. According to Mered- 
ith, 1 a 1% improvement in the maximum lift coefficient 
(or lift-to-drag ratio) could translate into an increased 
payload of 14 to 22 passengers on a large twin-engine 
transport. An optimal aerodynamic design of a multi- 
airfoil high-lift configuration requires careful consider- 
ation of both inviscid and viscous how phenomena. In 
particular, laminar- to- turbulent transition is a crucial 
issue for ground-to-hight scaling of high-lift how fields. 

Either of the familiar transition mechanisms for a 
single-element configuration are also relevant to tran- 
sition over a multi-airfoil configuration: those asso- 
ciated with streamwise instabilities in the form of 
Tolimien-Schlichting or Rayleigh modes, cross-flow 
and attachment line instabilities, and leading-edge 
contamination — a form of subcritical ( i.e ., bypass) 
transition. However, a unique transition mecha- 
nism in the case of multi-element how-fields involves 
the bound ary- layer contamination due to unsteady 
wake(s) from the preceding element (s) and/or addi- 
tional forms of vortical disturbances originating from 
the separated, cove flow underneath an upstream el- 
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ement (Fig. 1). While the “single- element” class of 
transition mechanisms have been widely studied in the 
literature, the wake-contamination issue has received 
little scrutiny thus far and is the focus of this paper. 

Interactions between turbulent w^akes and boundary 
layers have been the subject of much study. Most of 
the investigations, however, concentrated on the un- 
steady case that occurs, for instance, in turbomachin- 
ery, when the wake of an upstream blade impinges on 
a downstream blade; the unsteadiness of the impinge- 
ment region plays an important role in the dynamics 
of the how. Fewer studies can be found of the steady 
case. Squire 2 summarizes much of the investigations 
conducted prior to 1989. Particularly important is the 
work by Zhou and Squire 3 who examined the interac- 
tion between the wake of an airfoil and a flat plate. 
They found that, a region in which the wake and 
boundary layer are separated by a potential core is fol- 
lowed by a merging zone, in which the velocity profile 
in the outer layer is substantially different from that 
in a regular hat-plat a boundary layer. They also ob- 
served that the position of zero Reynolds shear stress 



Fig. 1 Sketch of the wake/boundary-layer inter- 
actions on an airfoil. 
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and that of zero mean-velocity gradient do not coin- 
cide, an important issue for eddy-viscosity turbulent 
models. 

Recently, Kyriakides and co-workers 4 '" 6 performed 
a series of experiments involving the interaction be- 
tween the wake of a large-aspect-ratio circular cylinder 
and a fiat-plate boundary layer. These experiments 
spanned a range of Reynolds numbers and cylinder 
diameters, including cases in which the wake still con- 
tained coherent eddies in the interaction region (i.e., 
the so called case of “strong interaction 55 ) and oth- 
ers in which the interaction took place sufficiently far 
downstream of the cylinder, such that the spanwise 
rollers had decayed significantly (i.e., “weak inter ac- 
tion 55 cases) . While they provided measurements for a 
range of Reynolds numbers, Kyriakides et ai. 4 " 6 did 
not clarify' how transition could occur (and turbulence 
could be self- sustained) at the lowest Reynolds num- 
bers examined. 

Due to the tremendous challenges inherent in a high- 
fidelity simulation of the wake-contamination problem, 
we also choose to study a simple, building-block prob- 
lem, namely that of interaction between a flat-plate 
boundary layer and the wake of a circular cylinder 
that is placed above the fiat plate. This model prob- 
lem provides a reasonable balance between simplicity 
(two-dimensional geometry involving a combination of 
two canonical flow fields) and the numerical challenges 
of a complex-geometry simulation including the wake 
generator and a boundary layer. Simultaneously, it al- 
lows one to investigate a range of issues, such as the 
effect of Reynolds number and the role of the spac- 
ing between the cylinder and plate (which allows one 
to traverse the continuous range from strong wake 
interactions to weak interactions that are analogous 
to conventional free-stream turbulence). Effects of 
flow' three-dimensionality and/or transverse pressure 
gradients (which lead to curved wakes) may also be 
captured, if necessary, by modifying the free-stream 
boundary conditions. Most significantly, a limited 
amount of experimental data is available for a range 
of flow and geometry parameters. 4 " 6 

In the present paper we report the initial re- 
sults from a series of numerical experiments de- 
signed to study the physical phenomena underlying 
wake/ boundary-layer interactions of this type. The 
first stage of this investigation involves the study of 
the strong interaction case, with the specific aims to 
clarify some issues left unresolved by the experimen- 
tal studies, and to document the flow field in detail. 
In particular, we will focus on the issue of transi- 
tion mechanisms at low' Reynolds number, w'hich was 
only partially addressed in the experiment. Thus, we 
choose to examine low Reynolds-number cases first, 
before attempting the more challenging simulations at 
higher Reynolds numbers. Future work will focus on 
the weak interaction case, as w f ell as other shapes of 


wake generators. 

Herein, we analyze direct and large-eddy simulations 
of two cases studied by Kyriakides et al , 4 having the 
same geometry but different Reynolds numbers, which 
lead to different responses of the boundary layer to the 
unsteady wake in the free stream. Despite the simplic- 
ity of the geometrical configuration, conducting DNS 
and LES of the above problem is a challenge for most 
numerical methods available today. In the present 
study we model the effect of the cylinder using an “im- 
mersed boundary” formulation. 11 This approach en- 
ables one to use codes in Cartesian coordinates, which 
present significant advantages, in terms of speed, accu- 
racy and flexibility, over codes that employ body'-fitted 
grids. 

The outline of the paper is as follows: the numeri- 
cal formulation of the problem including the geometry 
of the cylinder-plate configuration are described first. 
Then, the numerical methodology' will be presented. 
We will then present validation of the numerical tool 
in the context of tw'o building-block problems rele- 
vant to the flow configuration of interest; results of the 
simulations will be compared with numerical and ex- 
perimental data. Results for the wake/boundary-layer 
interaction are discussed next. Finally, conclusions 
and recommendations for future work are presented. 


Problem formulation 

In this work we present both DNS and LES re- 
sults. In the DNS case, the Navier- Stokes equations 
are solved: 


duj 

dxj 


-0, 


(i) 


duid 1 dp 

dt ' dxj UjUl ' pdxi 


+ i/V 2 Ui-r/i.(2) 


For the LES, we use the filtered equations of conser- 
vation of mass and momentum 


duj 

dxj 


-0, 


duj d 
dt dxj 




drji 

dxj 


J i 


( 3 ) 

( 4 ) 


(w'here the over-bar denotes filtered variables and the 
effect of the subgrid scales appears through "he SGS 
stresses — v~v~ - u{uj)\ f; and f i are the body 
forces used in the immersed-boundary method to en- 
force the no-slip conditions on the body. 

For both DNS and LES, the equations of motion 
axe solved numerically using a second-order accurate 
finite-difference method on a staggered grid. In all 
cases the grid is uniform in the spanwise direction 
y , stretched in the streamwise and normal ones (x 
and z, respectively) to allow accurate resolution of 
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Fig. 2 Sketch of the computational configuration 
and boundary conditions. The drawing is not to 
scale. 

the boundary layer and of the shear layers emanat- 
ing from the cylinder. The discretized equations are 
integrated in time using an explicit fractional time- 
step method, 7,3 in which all terms are advanced in 
time using the Adams-Bashforth method. The Pois- 
son equation is then solved using a direct solver, and 
the velocity is corrected to make the field solenoidal. 
For the LES, the SGS stresses are parameterized using 
the Lagrangian Dynamic eddy- viscosity model, 9 which 
has been shown to give accurate results in transitional 
flows. 

The effect of the cylinder, which does not coincide 
with the grid lines, is introduced using an “immersed 
boundary” formulation. In our implementation the 
enforcement of non-slip boundary conditions on the 
cylinder surface is performed using a method similar 
to the the “discrete forcing” approach introduced by 
Mohd-Yusof 10 and Fadlun et al. 11 With this approach, 
the solution is reconstructed at the grid nodes near the 
boundary to ensure satisfaction of the no-slip condi- 
tions. While Mohd-Yusof 10 and Fadlun et al. 11 used a 
simple one-dimensional reconstruction scheme, in the 
present study we employ a multidimensional proce- 
dure proposed by Balaras, 13 in which the solution is 
reconstructed along the line normal to the interface. 

The code is parallelized using MPI message-passing 
routines; the parallelization is performed by domain- 
decomposition process. Except for the Poisson solver, 
the computational domain is divided into a number 
subdomains in the z directions, each of which is as- 
signed to one processor. In the Poisson solver, after 
the velocity divergence is Fourier-transformed in the 
3 panwise directions, each processor is responsible :cr a 
set of span wise modes. 

We applied the following boundary conditions (see 
Fig. 2): 

1. At the inlet, a uniform velocity in the x direction 
w*as assigned. 

2. At the outlet, radiative boundary conditions 17 
were given. Furthermore, over the last 30% of 
the domain (which contains only 5% of the grid 


points) the mesh was stretched very significantly, 
to generate a buffer layer in which the equations 
are effectively parabolized because of the large as- 
pect ratio of the cells. 

3. In the spanwise direction, y , periodic conditions 
were used. 

4. In the freestream, we used a slip-wall boundary 
condition. 

5. On the cylinder surface, no-slip conditions were 
enforced using the methodology described in 
Balaras. 13 The cylinder was placed at a height of 
3.51? above the nominal leading-edge of the plate. 

6. On the z = 0 surface, w'e used ffee-slip condi- 
tions (w = 0, dujdz — dv/dz = 0) ahead of 
the plate, and no-slip conditions (ti* = 0) on the 
plate. A sharp transition between free-slip and 
no-slip was found to generate numerical oscilla- 
tions; therefore, a hyperbolic tangent profile was 
used to merge free-slip and no-slip conditions in 
a smooth fashion. As a result, the origin of the 
boundary layer on the plate cannot be determined 
precisely; in the present calculations, the virtual 
origin of the boundary layer appears to be ap- 
proximately 1.57? upstream of the actual leading 
edge. 

Two calculations were carried out, corresponding to 
cases 1 and 2 of Kyriakides et al. 4 The geometry for 
the two cases was the same, but the free-stream ve- 
locity was varied to obtain Reynolds numbers Reo 
equal to 385 and 1,155, respectively. The definition 
of the Reynolds number is based on the free-stream 
velocity and cylinder diameter. The extent of the com- 
putational domain in the streamwise and wall-normal 
directions was 1001? and 301?, as shown in Fig. 1. The 
spanwise width of the domain was 2ttD wide, a length 
sufficient to include several rib vortices, and to resolve 
the three-dimensional structures in the cylinder wake. 

Code validation 

The code used in this paper has been exten- 
sively validated for a variety of turbulent 14,10 and 
re-laminarizing 16 flows. The “immersed boundary” 
formulation which is used to enforce non-slip boundary 
conditions on the cylinders surface, has also been ex- 
tensively tested in the context of the present code. 13 In 
particular, a series of three-dimensional computations 
were conducted of the flow around a circular cylin- 
der at Reo — 300. These involved a total number of 
grid points ranging from 1.5 x 10 6 to 4.7 x 10 6 . The 
agreement of the results with reference computations 
using boundary fitted coordinates demonstrated the 
accuracy of the method. Details on the above compu- 
tations can be found in Balaras. 13 
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Fig. 3 Streamwise development of the skin-friction 
coefficient Cf. 



Fig. 4 Profiles of the streamwise turbulence in- 
tensity, Urms- 

The grid required for accurate simulations of bypass 
transition needs to be determined carefully. To this 
end, we performed simulations of bypass transition on 
a flat plate due to free-stream turbulence. We used the 
setup adopted by Yoke and Yang. 18 The calculation 
was started at the streamwise location Re x ~ 6, 000; 
the inflow’ condition consisted of a Blasius boundary 
layer, on which turbulent fluctuations were added. 
The fluctuations were obtained from a separate calcu- 
lation of spatially developing homogeneous isotropic 
turbulence. To avoid introducing excessively large 
fluctuations near the wall, the turbulent fluctuations 
were windowed in such a way that they vanished at 
2S 0 /3 (where 5 0 is the inflow boundary-layer thick- 
ness). Two cases were run, with a 3% and a 6% 
freestream turbulence amplitude, to match the calcu- 
lations 18 and the experimental data. 19 The descrip- 
tion of the inflow conditions both for the experiment 
and the LES calculations was not, however, detailed 
enough to allow us to reproduce exactly the experi- 
mental setup. 

We performed DNS of both cases. For the 3% case, 
we 850 x 96 x 96 points to simulate a domain that 
had dimensions 714 S 0 x 32 S Q x 48d 0 ; For the 6% case. 
550 points were used in the streamwise direction to 
discretize a box that was 476 5 0 long. The mesh was 


uniform in the spanwdse direction, and stretched in the 
x and r directions. In the transition region, the grid 
size was approximately Ax = 0.2 8, A y — 0.08(5; at 
least 40 points were used to discretize the local bound- 
ary layer thickness, <5. For the LES of the 6% case, 
300 x 64 x 96 points were used, resulting in coarser 
resolution in x and y: Ax = 0.4<5, A y — 0.1 5. 

The streamwrise development of the skin-friction co- 
efficient Cf — 2r w jpUl Q IS shown in Fig. 3. One can 
observe a higher skin-friction coefficient in the lami- 
nar region, due to the perturbation introduced at the 
inflow, which generates a pressure disturbance that 
propagates w~ell inside the boundary layer and dis- 
torts the mean- velocity profile in the near- wall region. 
It was verified that, in the absence of inflow distur- 
bances, the Blasius profile was recovered. The slope 
of the skin-friction coefficient distribution is, however 
correct. The transition begins at the correct location 
for both values of the freestream turbulence intensity, 
and the extent of the transition region is also predicted 
7v*ith reasonable accuracy. The LES calculation pre- 
dicts the a slightly delayed beginning of the rise of 
Cf , and a somewhat more rapid transition; the results 
are, however, in reasonable agreement wdth the DNS 
and the experiments. In Fig. 4 the distributions or the 
r.m.s. turbulence intensities at Re x = 290, 000 (for the 
3% freestream turbulence case) are compared with the 
experimental data of Roach and Brierley. 19 The agree- 
ment is again satisfactory, given the uncertainty in the 
inflow condition generation. 

Finally, a grid refinement study was carried out on 
the real geometry. Two grids were tested, a coarse one 
that used 576 x 48 x 192 grid points and a medium one 
with 864 x 72 x 288 points. Calculations on a finer grid 
with 1056 x 128 x 384 points are currently underway. 
The coarse grid was similar to the one that gave grid- 
converged results in the cylinder calculation. Figure 5 
shows the mean velocity profile and the trace of the 
Reynolds stresses q 2 = (u^u') (throughout this paper 
(■) represents a long-time average, and /' = / - (/) 
is the fluctuating component of /) at three locations 
dowmstream of the cylinder. Very little difference cam 
be observed between the coarse and medium grids for 
the mean velocity profiles; the coarse resolution of the 
wake and boundary’ layer results in some difference in 
q 2 between the coarse and medium meshes, but the 
agreement is altogether satisfactory. In the follow’ing, 
the medium grid results will be presented. 

Notice that the resolution of the boundary layer in 
the medium grid compares favorably with the bypass 
transition case mentioned above. In the transition re- 
gion, the grid spacing was Ax = 0.07<5, Ay — 0.06(5 
(where d is the local boundary layer thickness). At 
least 80 points were used to discretize the boundary- 
layer thickness. 
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Fig. 6 Mean velocity contours. 


Rcd = 385; Ren = 1, 155. 
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Fig. 5 (a) Mean velocity profiles and (b) profiles 

of q 2 , Rcd — 385. 576 x 48 x 192 grid; 864 x 

72 x 288; x Blasius solution. 
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Fig. 7 Streamwise development of: (a) the ve- 

locity at the edge of the boundary layer; (b) the 
acceleration parameter X. 

Results and discussion 

The mean-velocity contours from the 
wake/boundary-layer simulations at Reo = 385 

and Re o — 1, 155 are shown in Fig. 6. As expected, 
both the wake and the boundary layer are thinner 
at the higher Reynolds number. No merging is 
observed within the computational domain. Due to 
the presence of the cylinder, the boundary layer near 
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Fig. 8 Stream wise development of statistical quan- 
tities. (a) skin friction coefficient C/; (b) Reynolds 
number based on the momentum thickness, Re#; (c) 
peak 3 1 reamwise and spanwise Reynolds stresses in 
the boundary layer. 

the leading edge first encounters an accelerating outer 
flow, and then a decelerating one. The corresponding 
acceleration parameter K — {yj ) (dU^ fdx) (where 
Ux> is the velocity at the edge of the boundary layer) 
is initially large and positive. The region of large 
positive values of K is, however, limited to the region 
x/D < 1 (Fig. 7); over most of the Hat plate the 
acceleration parameter is below the value that may 
result in re~!aminarization (K ~ 3 x 10~ 6 ). In fact, 
for x/D > 22 the boundary layer is subjected to a 
mild adverse pressure gradient that is expected to 
favor transition to turbulence. 

The streamwise development of several statistical 
quantities is shown in Fig. 8. The skin friction coeffi- 
cient Cf (Fig. 8a) from both calculations is very close 
to the value corresponding to the Blasius boundary 
layer. In the high Reynolds-number case we observe a 
significant increase of Cf (towards the end of the com- 
putational domain); this issue will be discussed further 
below. The Reynolds number based on the momentum 
thickness, Reg, is shown in Fig. 8b. For both calcu- 
lations Reg remains very low. As mentioned before, 
the Reg values achieved in the Re& = 385 case are 
unquestionably below the minimum values at which 
self-sustained turbulence is typically observed in flat- 



Strouhal 


Fig. 9 Power spectra of the streamwise velocity 
fluctuations near the wall, (a) Data of Kyriakides 
et a/.; 4 (b) present calculation, Rcd — 385. 

plate boundary layers. The Reg values from the high 
Reynolds-number calculation towards the end of the 
computational domain are perhaps marginal in this 
regard. 

The power spectra of the streamwise velocity at 
z/D = 0.15 are compared in Fig. 9 with the exper- 
imental data of Kyriakides et o/. 4 The experimental 
data was taken from Fig. 7 of the paper cited, which 
did not specify the distance from the wall at which 
the data was measured, nor the scale of the ordinate; 
the frequency has been rescaled to yield a dimension- 
less Strouhal number. Our calculations are in general 
agreement with the experimental observations; the 
shedding frequency is 0.22 (in the experiment a value 
of 0.21 was observed) and we observe a filling-up of the 
spectrum at both lower and higher frequencies with 
distance downstream of the cylinder. We do not ob- 
serve, however, the harmonic peak measured in the 
experiment at Si = 0.42. It is unclear whether this 
discrepancy is due co experimental or numerical er- 
rors, to insufficient sample convergence in the DNS, or 
to a difference in the location of the measurements. 

Kyriakides et al A observed that time-traces of the 
velocity at this height had a sinusoidal behavior up 
to x/D ~ 7.3 and 2.7 for the low and high Reynolds- 
number, respectively; they identified transition with 
the loss of the sinusoidal behavior. The time histories 
of the velocity fluctuations (not shown) show depar- 
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Fig. 10 Mean velocity profiles at various locations. 

Reo ~ 385; Ren = 1, 155; x B Iasi us profile. 

Each profile is shifted by 1 unit for clarity. 

ture from the sinusoidal behavior at similar locations. 
We believe, however, that a quantitative definition 
of the onset of transition is more desirable than the 
subjective one adopted by Kyriakides et ah 4 One pos- 
sibility is to define the location of the transition point 
in terms of some threshold value of the turbulent fluc- 
tuations. Figure 8c shows the developments of the 
peak (u'u') and (vV) Reynolds stresses, normalized 
by the local wall stress r w . The streamwise Reynolds 
stresses reach values comparable to those expected in 
turbulent flows in the Reo = 1, 155 case, while they 
remain at about half of this level when Reo — 385. 

An important measure of transition is the onset of 
three-dimensionality. We observe that the spanwise 
Reynolds stress (v'v') begins growing earlier for Reo = 
385 than for Reo — 1, 155. This is a consequence of 
the faster spreading of the cylinder wake at the low 
Reynolds-number. It wall be showm later that the first 
onset of three-dimensionality in the boundary layer is 
caused by the passage of the vortices shed in the wake 
of the cylinder. The onset of three-dimensionality in 
the boundary layer is very well correlated with the 
onset of three-dimensionality in the wake. 

Figures 10 and 11 show the profiles of the mean ve- 
locity and q 2 at several locations downstream of the 
cylinder. We observe a clear separation between the 
boundary layer and the wake near the leading edge, 
which disappears as one moves downstream. Up to 
x/D = 20 the velocity profiles agree very well with 
the Blasius solution at both Reynolds numbers. In the 
low Reynolds- number case wrn observe a mild growth 
of fluctuations in the boundary layer (evidenced both 
by the streamwise development of (u'u'), Fig. 8, and 
by the profiles of q 2 , Fig. 11). These fluctuations, how- 
ever, appear to be due entirely to the advection caused 
by the vortices shed by the cylinder, wdiich, as they are 



Fig. IX Profiles of q 2 at various locations. 

Rea — 385; Reo = 1, 155. Each profile is shifted 

by 0.05 units for clarity. 



Fig. 12 Profiles of q 2 in wall units at various lo- 
cations. Reo — 385; Reo = 1,155. Each 

profile is shifted by 5 units for clarity. 


convected downstream, induce ejections of low- speed 
fluid from the near-w-all region. These velocity fluctu- 
ations, in fact, do not result in the levels of Reynolds 
shear stress ( u'w ') typical of the turbulent fiat-plate 
boundary layer (see belowr). 

The high-Reynolds-number calculation shows a dif- 
ferent behavior. Significant differences appear in the 
velocity profile at the downstream locations. At the 
last location showm ( x/D — 55), Re x ~ 58,000 and 
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Fig. 13 Profiles of (uw f ) in wall units at various 

locations. Reo = 385; Ren = 1, 155. Each 

profile is shifted by 1 unit for clarity. 

Res * ~ 340. Although the flow has not achieved a fully 
developed turbulent state yet, there are several indica- 
tions that transition to turbulence is taking place. The 
turbulent kinetic energy (Fig. 11), for instance, shows 
a near- wall peak that is increasing and moving towards 
the wall. The value of this peak and its location in 
wall units are those expected in a turbulent flat-plate 
boundary layer (Fig. 12). The Reynolds shear stress 
{u'w*) (Fig. 13) exhibits a similar behavior: its peak 
value and distribution for the high-Reynolds-number 
case are similar to those expected in turbulent flat- 
plate flows, while in the Reo ~ 385 case the peak is 
significantly lower. 

The Reynolds stresses appear to move towards a 
turbulent state more rapidly than the mean velocity 
profile; the shape factor H — 5* / 9 remains close to the 
laminar value of 2.6 throughout the computational do- 
main, and the mean velocity profile has not achieved a 
logarithmic profile by the end of the computational do- 
main (Fig. 14). However, significant levels of turbulent 
fluctuations (r.m.s. levels between 1 and 10% of the 
freestream velocity) are observed within the bound- 
ary layer. At the downstream stations, the ffeestream 
turbulence at the edge of the boundary layer is very 
significant, between 3 and 5% of the velocity at the 
boundary-layer edge (for the low and high-Reynolds 
number cases, respectively), levels that are able to 
trigger bypass transition in boundary layers. Com- 
putations that use a longer domain should be carried 
out to determine the exact route to turbulence in the 
present configuration. 



Fig. 14 Mean velocity profiles in wall units at 

various locations. Reo — 385; Reo = 1, 155. 

Each profile is shifted by 10 units for clarity. 

Figures 15 and 16 show instantaneous contours of 
the streamwise velocity fluctuations, and iso-surfaces 
of the second invariant of the velocity gradient tensor, 

(5) 

where Sij is the strain-rate and f2* ; the rotation-rate 
tensor. In regions where Q > 0, the vorticity is due to 
rotational motions, rather than to shear. 

In the low-Reynolds-number case (Fig. 15) the cylin- 
der wake is very w’ell-defined; staggered row’s of quasi- 
2D vortices are shed, which are joined by clearly ob- 
servable rib vortices. Three-dimensionalities quickly 
develop. There is a high correlation between the ap- 
pearance of three-dimensionalities in the rollers, and 
the appearance of low- and high-speed streaks in the 
boundary layer (a very clear example can be observed 
around x/D — 20). This feature is connected to the 
increase of the spanwise Reynolds stress (uV) men- 
tioned above. These streaks are convected with the 
freestream velocity, confirming the observation that 
they are induced by the outer flow” in turbulent flows, 
the typical near- wall streaks, which are generated in 
the inner layer, are convected at speeds that are 60- 
70% of the freestream velocity. 

Another effect of the shed vortices on the near-wall 
region can be observed in Fig. 15c. Corresponding 
to the vortex street, one can observe alternating re- 
gions of high- and low’-speed fluctuations, due to the 
advection from the upper and lower row’s of vortices, 
respectively. These regions correspond to Q2 and Q4 
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= 0.2 plane. Reo = 385. (a) Prospective view; (b) top view; (c) side view. 
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events, 20 and are responsible for regions of significant 
correlation u'w 1 . As the rollers become weaker, no new 
streaks are generated, and the flow does not appear to 
evolve into a fully turbulent state, the main distur- 
bances being due to the streaks that are convected 
downstream, and weaken progressively. 

At the high Reynolds- number (Fig. 16) a some- 
what different picture emerges. Smaller-scale vortices 
can be observed in the wake, consistent with the in- 
creased Reynolds number. Also, the wake breaks down 
more rapidly. Well-defined tollers and rib vortices can, 
nonetheless, be observed. Some low-speed regions in 
the boundary’ layer can still be observed, correspond- 
ing to the onset of three-dimensionality in the rollers; 
unlike the low-Reynolds-number case, however, these 
streaks do not weaken significantly as they are con- 
vected downstream. The decreased viscous dissipation 
due to the increased Reynolds number certainly plays a 
role in this, but a self-sustaining turbulence generation 
cycle may also be beginning to appear. This conjecture 
is supported by the appearance of quasi-streamwise 
vortices in the near-wall region for xjD > 50- 

Conclusions 

The work presented in this paper had a two-fold 
objective, namely, to establish the accuracy of the 
numerical code for wake-boundary-layer interaction 
problems and to clarify unanswered questions regard- 
ing prior experiments 4 related to such interactions. We 
believe that both of these goals have been nearly met 
through the work described here. In addition to code 
validation for building-block problems involving the 
isolated problems of cylinder wake and bypass tran- 
sition in a fiat-plate boundary layer, we have carried 
out large-eddy and direct simulations to understand 
the details of the strong interaction between a cylinder 
wake and a fiat-plate boundary layer at two Reynolds 
numbers [Rep ~ 385 and 1,155) selected from the low- 
Reo range in the experiment, for which transition was 
apparently observed. 

We observe a region in which the cylinder wake 
perturbs the boundary layer by advecting low-speed 
fluid from the wall into the outer region, and high- 
speed fluid from the core towards the wall, as rollers 
of alternating rotation are convected above the bound- 
ary layer. Three-dimensional disturbances appear in 
the boundary layer in correspondence to kinks in the 
spanwise rollers. The appearance of stream wise vor- 
ticity in the wake is visually well-correlated with the 
presence of low- and high-speed streaks near the wall. 
This, and the fact that these streaks are convected 
downstream at the freestream velocity (rather than a 
convection velocity typical of near-wall phenomena) 
indicate that the perturbations in the boundary layer 
are forced from the outside, and not generated within 
the boundary layer itself. As the wake decays, the 
advection effect becomes less significant, but the wake 


generates high levels of freestream turbulence, ~hat are 
expected to affect the How downstream. In the low- 
Reynolds-number case the flow tends towards a more 
quiescent state towards the end of the computational 
domain and the statistics never reach turbulent levels. 
In the high-Reynolds-number case, on the other hand, 
the appearance of coherent, quasi-streamwise vortices 
is observed, which might be the prelude to the es- 
tablishment of a self-generating cycle of sweeps and 
ejections. This was supported by the establishmnent 
of Reynolds-stress levels comparable to those observed 
in turbulent flows. 

While seme features of the computed flow fields are 
in satisfactory agreement with the measured data, the 
simulations indicate that transition does not really 
take place at the lower Reynolds number: the bound- 
ary layer is simply buffeted by the unsteady Karman 
vortex street shed by the cylinder. For Red = 385 the 
fluctuations generated in this manner cannot be self- 
sustained due to low-Reynolds-number effects, and the 
boundary-layer flow does not transition to turbulence. 
We believe that the cause behind the discrepancy 
between the observations of Kyriakides et a/. 4 and 
the above results is to be found in an unsatisfactory 
criterion for transition detection in the experiment, 
and in the lack of sufficiently detailed measurements 
of the type particularly required for transition stud- 
ies. In contrast, in the higher Reynolds-number case, 
boundary-layer fluctuations persist even after the wake 
has decayed. This can be attributed, in part, to the 
spatial growth of the boundary layer during the in- 
teraction, which brings the Reg values closer to the 
typical range of minimum critical Reynolds numbers 
required for self-sustained turbulence in zero-pressure- 
gradient flows; some evidence could be observed that 
a self-sustaining turbulence generation cycle was be- 
ginning to be established. 

Near-term extensions to this work will involve addi- 
tional simulations to clarify the flow behavior down- 
stream of the present domain. We also plan to inves- 
tigate the possibility of higher- He simulations, which 
will also allow the present study of strong wake inter- 
actions to be extended to weak interactions. It will 
be interesting to find out the extent of similarities be- 
tween the weak interactions and the familiar case of 
bypass transition due to moderate levels of free-stream 
turbulence. 
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